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Abstract. We investigate the third order spontaneous parametric down- 
conversion process in a nonlinear media with inversion centers. Specifically, we 
analyze in details the three-photon differential count rate in unit frequency and 
angular regions, total count rate and measurement time for rutile and calcite 
crystals which have comparatively large cubic susceptibilities. Special attention 
is given to consideration of limited frequency and angular detection ranges in 
order to calculate experimentally available detection rate values. 
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1. Introduction 


Generation of photon-number (Fock) states of light 
is one of the main tasks in quantum optics. They 
are interesting not only from fundamental, but also 
from practical points of view because of their necessity 
for solving problems of quantum communications and 
linear optical quantum computations. While problems 
of single-photon and biphoton-state generation are well 
studied, the direct and non post-selective generation of 
higher-order Fock states is still an attractive challenge. 

In this work we consider the problem of three- 
photon state generation. Non-classical properties 
of such states enable heralded emission of photon 
pairs [TJI^ as well as preparation of three-body 
entangled states (for example, Greenberger-Horne- 
Zeilinger (GHZ) states mil])- 

There are several proposed solutions for the 
problem of three-photon generation such as cascaded 
or postselective second-order nonlinear processes [6 ■ 12 
and formation of approximate photon triplets by SPDG 
photon pairs together with an attenuated coherent 
state [^. All these approaches give relatively low 
photon generation rates (up to 45/minute [^) and have 
a big contribution of low-photon-number impurities. 

On the other hand, the most natural way 
to generate three-photon states is a third-order 
spontaneous parametric down-conversion (TOSPDG). 
Unlike the other techniques, it enables to generate a 
three-particle entanglement in continuous degrees of 
freedom, such as energy and momentum. This problem 
was previously studied theoretically, but to the best of 
our knowledge no experimental results were reported 
for direct spontaneous generation of triplets based on 
Only stimulated third-order parametric down- 
conversion was demonstrated by seeding triplet modes 

(Hum 


There are two approaches for TOSPDG genera¬ 
tion: in bulk crystals [6l [l6}|l8] or in optical fibers 
19 -24 . The bulk crystals allow satisfying simply the 


phase-matching condition using different polarization 
modes, but a spatial multimode structure of three- 
photon light and a limited crystal length complicate 
a high photon conversion probability and an effective 
detection. One can increase the interaction length and 
decrease the spatial mode number by using optical 
fibers. In this case phase-matching condition can be 
realized while the pump and three-photon light prop¬ 
agate in different spatial modes or by using the quasi- 
phase-matching. But a small mode overlap (~ 10“^ 
and a high absorption coefficient for the pump 
(in a visible and especially in UV case) also limit the 
generation rate. 

Our work presents a theoretical description of 
the third-order parametric down-conversion in crystals. 
Special attention will be paid to crystals with inversion 


centers, having zero and hence prohibiting all 
three-wave processes. That is of great importance 
especially in the case of triple generation with seeding 
beams because three-wave processes are much more 
intense and may suppress generation of triplets as well 
as their detection. 

The paper is organized as follows. In Sec. 2 we 
evaluate the three-photon count rate in unit ranges 
of frequency and transverse wave vector and the 
integral count rate over all detectable frequencies 
and transverse wave vectors of scattered photons in 
a collinear degenerate regime of generation for the 
type-I and type-II phase-matching. In Sec. 3 we get 
estimates of the minimal measurement time sufficient 
for distinguishing signal triple coincidences from noise 
ones. Then in Sec. 4 our estimates are specified for 
two nonlinear crystals with inversion centers: calcite 
and rutile. And finally in Sec. 5 we discuss obtained 
results. 


2. Calculation of photon count rate 


As the process of third-order SPDG (TOSPDG) is 
similar to two-photon SPDG, in this section we follow 
the approach developed by D. N. Klyshko for biphotons 
in [^, though somewhat extended for the case of 
triplets. 

Let a pump photon in the mode kp, uip be decaying 
for three photons in modes fci,wi, fc 2 ,W 2 and ks^Lo^ 
and let a pump be a monochromatic plane-wave 
propagating along the z-axis. Under these assumptions 
the photon energy and transverse momentum are 
conserved and the non-conservation of the longitudinal 
momentum determines the phase mismatch Afc^: 


Wi -|- UJ2 + 0J3 — OJp — 0, (1) 

gi + <72 + <73 = 0, (2) 

kiz + k2z + hz - kp = Ak^, (3) 

where denote perpendicular components of ki. 

In the second order of the perturbation theory 
TOSPDG is described by the Hamiltonian 


H = ^ f d^r X^^'’ck,Ck,Ck^aiy^^al^Ep x 


V fcl,fe 2 .fe 3 


exp 


i{kp — ki — k 2 — k^)?— i{ujp — wi — W2 — ^3)^ 

+ H.C., (4) 


where Ep is the amplitude of the pump considered 
as a classical monochromatic plain wave, V is the 
interaction volume. 




2'KfUjJk 


( 5 ) 


V is the quantization volume and aj(,j are the photon 
creation operators for modes ki (i = 1, 2, 3, p). 
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The rate of transitions per unit spectral and 
transverse-wave-vector ranges of each photon diVi and 
dtfi is determined by the Fermi Golden Rule 

/? - - - — r PW sinc^ (X 

^•■u;igiti;2(32‘^393 ^ I ^ j 


Here . .]o means that the expression in brackets is 
evaluated at the exact collinear degenerate regime of 
generation. All mixed derivatives are equal zero. 

Let’s consider separately two types of phase¬ 
matching, type-I and type-II. 


, ( 6 ) 


where 

= S{qia: ■ 
— aj 1 ajp /3 


T = h 


(3) 


q2x + q3x)S{qiy -\- q2y + q3y), 

(5(r2) = 5(r2i - 1 - H 2 -f H 3 ), 

LUiUJ2Ui3f (c^{2Tr)‘^nin2n3np) , (7) 

is the pump power 


ni are the refractive indices, Wp 
and c is the speed of light. 

To calculate the total three-photon generation 
rate, we have to integrate the differential rate of Eq. ([^ 
over spectral and angular regions, restricted by features 
of the detection scheme: 


Rt « ITl Wp, where 
I = sine 


( 8 ) 


2 1 6^^Hq)S{n) X 


X duji du}2 duj3 dqi dq2 dq^ (9) 

The above-mentioned restrictions of the integration 
ranges will be applied to variables of all three 
photons. This is the approach needed for description of 
experimentally measurable correlations of photons. In 
contrast, in earlier publications [T6p7 filtering related 
to the features of detectors was assumed to be applied 
only to variables of one photon of a triplet whereas 
for two other photons ranges of integration were taken 
unlimited, which is sufficient only for description of the 
single-photon count rate. 

In the biphoton case one can calculate the total 
count rate, integrating along the phase-matching 
curve, defined by the equation A fcz( 91,(72 = 

—( 7 i,Hi,H 2 = —Hi) = 0. This case is comparatively 
simple as there are only two independent integration 
variables. Below we will derive the expressions for 
the curves analogous to the biphoton phase-matching 
curve in the case of triplets for the type-I and type-II 
TOSPDC collinear degenerate regime of generation. 

In order to do this, let us expand Ak^ in powers 
of Hi, qix and qiy up to the second order: 

kzi = \Jk1 - qj, ki = n{LOi)^, i = l, 2 ,3, ( 10 ) 

Akz = kzi + kz2 + kz3 - kp = [ki -|- -I- ^ 3 ]^ - kp-\- 


E 

z=1.2,3 


dkj 

flqix 


qix 


dkj 

dqiy 


qiy 


dkj 

duji 


Hi 


1 

Id'^kil 

2 , 1 


2 , 1 

to 

2 

do^ 

L J 

9jx + 

0 ^ 

1 

Ty + 2 

0 



( 11 ) 


2.1. Type-I phasematching 


In the case of phase matching of the type-I (e—J-ooo) 
the first-order derivatives of ki in Eq. ( [To| appear to 
be identical for different i and, as ^ qix,y = 0 and 




0, Eq. (Ill takes the form 


I 


Hs) = ^H| - aq^, 


( 12 ) 


where 


a = -- 


d'^Akz 

~W 


-h 

- ^Kp, 


/3 = 


9 I = 


1 

'd^Ak^' 

1 


2 

[ deaf \ 

Q Jttc^ 

[ 3A?J 


2 = 1 , 2 ,3 


ni 


= E nl 

i=l,2,3 


J 0 


Thus, Akz depends on two variables, q^: and Hs, 
and hence, the exact-phase-matching curve is defined 
by the equation /3H|. — aq§^ = 0 (the dotted green line 
in Fig. la). 

Now, let us change the integration variables in 
Eq. ^ Hi, H 2 , H 3 Ha, Hb, He (Fig. Ic): 

2 H 3 — Hi — H 2 Hi -|- H 2 

cs ■ 

Hi -|- H 2 -l- H 3 

= — 75 — 

and 9 i, 92 , 93 ^ <74, qs, qc- 


293 -91-92 ^ 9 i + <72 

“ = — 75 —' 

^ 91 + 92+93 

In these variables i5(H) = 5(He) and 5^^\q) = 
5^‘^\qc) in owing to which integrals over He 
and qc are easily taken. As for the other variables, 
let us introduce the polar coordinates in the planes 
(qAx.qAy), {qBx,qBy), iqA,qB), and (Ha,Hb): 

QAx = QaCOSC^a, QAy = qASm(j)A, 

qBx = qB cos 4>b , qBy = qB sin (fs, 
qA = 9s cos 4>q, qB = 9s sin (j)q, 

Ha = Hs cos (/>n, Hb = H^ sinc/iQ. 

The Jacobian of transformation to polar coordi¬ 
nates is given by J = H-^q^ cos (j)q sin (j)q and Eq. (|^ 
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Figure 1. (a) Function sinc^ [Akz{qj:-, and corresponding spectral distributions in qA — qs parameters (b) and in Q.a — 

parameters (c) for type-I TOSPDC. Green dotted line in (a) corresponds to the phase-matching curve defined by the equation 
Akz{qY,i = 0- Calculation is made for 10 cm long rutile crystal at pump wavelength 532 nm. 


takes the form: 


7 = / dqY;d^-^d(j)qd(j)Adcj)Bd(j)n J 


Akzl 


{2nr 


sinc^ 


/AkJ 


(13) 


The phase-matching area in the plane is 

shown in white in Fig. la. Each point at the exact- 
phase-matching curve in this area corresponds to one 
quarter of a ring of the spectral distribution in qA — qs 
variables (Fig. lb), and to a full ring of the distribution 
in the — fig variables (Fig. lc).The integrals over 
the qA — qs and Qa — distributions are included 
into the Jacobian J. 

The integral over qs in Eq. (13 \ can be 
approximated by the product of the integrand at the 
exact-phase-matching curve with the width of the 
phase-matching area in the qe direction, Aq-^. The 
latter can be found from the equation : 

(/3n| - a{q^ -I- AqY.f)l/2 = tt ^ q^ Aq^ = 

la 


(14) 


which reduces Eq. (13) to the form 

2 al a a^l 


Umax 5 


(15) 


where 

^Emaa; = min 


27rc 


Up uup 


27rc knO 


p^max 


^min 3 3 ^max 3 

Here we have taken into account the limitation of the 
spectral (from Xmin to Xmax) and angular (no more 
than 9m,ax) detection range. 


2.2. Type-II phasematching 

Let us consider now the TOSPDC process with 
the type-II (e—^ooe) phase-matching. Let the optical 


axis of a crystal is located in the {xz) plane. In this 
case the first-order derivatives in Eq. are not equal 
and the decomposition of Akz has the following form 
(we assume that the photons 1 and 2 are ordinary and 
the photon 3 is extraordinary): 


Akz — + ^^ 2 ) + /de^3 + (^ed3x + 

^loxiqix + Q2x) + 2'^exq3x + 

-jloyidly + dly) + -jleyqty, (16) 

where 

O _ ^^1.2 a - dks _ dks 

Po — o 5 Pe — Q 5 OLe — ^ , 

0101^2 oqsx 

d'^ki 2 

We took into account here that inside a small angular 
range in the a;-direction the terms with the second- 
order derivatives 7 o,ex are much smaller than the terms 
with the first-order derivatives Q:o,ej but in the y- 
direction dAkz/dqiy=0 and, hence, the second-order 
derivatives 'yo,ey have to be retained. 

By denoting 7 ~ ~ ^oy, we get 


Akz — {Po — Pe){^l + D 2 ) — Q:e(9la; + q2x) + 

lidly+dly). (17) 

and the integral I of Eq. takes the form 




dVt+dQ.-dq+dq-, 


where 


zb fl2 qix i q2x 

a+ = y/2ae, P+= V2{Po - Pe), 


with the intervals of qf 2 y where | < 2tt estimated 

as Aql = 27r/Z|7|. 
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Figure 2. (a) The function sinc^ [Afcj (0+, qj^)l/2\ and the integration areas over the variables qi, q 2 (b) and over the variables f!i, ^2 
(c) for type-II TOSPDC. Green dotted line in (a) corresponds to the phase-matching curve defined by the equation Akz{fi+, q+) = 0. 
Calculation is made for a 1 mm long calcite crystal at the pump wavelength 532 nm. 


Hence, the exact-phase-matching curve for type- 
II TOSPDC has the form /3+ri+ — a+q+ = 0 (the 
green dotted line in Fig. [^). Each point of the phase¬ 
matching area { 0 +, g+j corresponds to the ranges A( 7 _ 
(Fig.[^) and AO_ (Fig.|^), which can be found from 
phase-matching conditions: 

Ag_ = 2{q^V2- | 9 +|) = 2 

(shown as blue arrow on Fig. 

AO_ = 2{nmV2- |0+|) 

(shown as red arrow on Fig. §)• 


3. Evaluation of the measurement time 

The measurement time T 3 can be defined as the time 
sufficient for extracting the signal triple coincidence 
count rate Rs from noise . Mathematically this 
means that the total number of triple coincidence 
counts (found as the difference of two measured rates, 
of the sum of signal and noise counts and, separately, 
of only noise counts) is at least tc,oo 

times bigger than its standard deviation CTs , where 
tc,oo is the Student’s t-factor for a confidence level C. 
Taking into account the Poisson distribution of photo 
counts one can calculate the dispersions: 



Similarly to the type-I case we approximate the 
integral over by the product of the integrand at 
the exact-phase-matching curve with the width of the 
phase-matching area Aq_(_ = At:/\ ae\l and get 




d^l+ ■ AD_Ag+A(7_ = 




J dn+■ {n^V2 - \n+\) X 



where Tlmax, and q^ are determined by detection 
ranges in frequency and angle shown in Fig. 2: 


Qm — 


Tim = min 


Tlm.ax = min[n^, 12 ,], with 


2ttc 


LUp UJp 

T’ T 


2ttc 


O = ^ 

“q — qm p ■ 


= y(2i?(f)+i?®)T3 

So, we obtain the equation for T 3 : 
iV(3) =T3i?i3) = tc,ooX = 

With given quantum efficiency rj and the noise 
count rate Rn'^ of each detector^ | (we assume 
that all detectors have approximately the same 
characteristics), the temporal resolution of electronics 
(typically limited by a detector jitter) 5t, and R^ 
evaluated in ([^, we get 

R^^^ = (19) 

i?(3) = (20) 

^ We consider an ideal case when the number of the noise counts 
equals to the number of intrinsic detector’s dark counts and 
the light noise is neglected. 
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where the parameter ^3 (and also the parameter ^2 “ 
see below) characterizes features of non-polarized beam 
splitters to be used in a possible experimental setup for 
dividing TOSPDC signal into three channels. Hence, 


2 


T3 = t 


2 

C,oo 


St^ + 


( 21 ) 


Similar expressions can be derived for the minimal time 
T 2 (Ti) required for distinguishing two-photon signal 
coincidence and single-photon counts from the noise: 

, 2 


T2=f 


C,oo^ 


{^Riy dT+^2RTv‘^ 


Ti =t- 


2R!'n^5T 


C,oo 


(?2 -Rt??^) 
RtV 


(Rtv) 


( 22 ) 


(23) 


and in case of two consistent 30/70 and 50/50 beam 
splitters, which provide approximately equal power in 
three channels, parameters ^ 3^2 are given by ^3 = 0.22 
and ^2 = 0.75. 


4. Example: calcite and rutile crystals 


Let us make estimates for two specific crystals, 
calcite and rutile, having comparatively large cubic 
susceptibilities. For these two crystals and for four 
pump wavelengths, 266, 325, 405 and 532 nm, the 
results of calculations are presented in Table B By 
using data about crystal refractive indices of Refs. 

27 , we found values of angles between the optical 
axes of crystals and the pump propagation direction 
providing collinear emission of TOSPDC photons. 
For these orientations of crystals, with the use of 
matrix elements determining 


28- 


31 


and with 

the dependence of on the angle between the pump 
propagation direction and the crystal optical axis |32| 
taken into account, we calculated values of the effective 
cubic susceptibility X^eff both crystals and for 
all collinear TOSPDC regimes indicated in Table 
Together with X^eff’ present in Table the total 
count rates® , calculated from Eq. for rutile and 
from Eq. (18 1 for calcite. 


Note that though rutile is a positive crystal and in 
the type-I phase-matching all three TOSPDC photons 
are not ordinary, it can be shown that even in this case 


the calculation based on (151 results in inaccuracy of 
Akz about 10 percents. 

The phase-matching conditions in calcite at the 
mentioned pump wavelengths are satisfied for each 

type of phase-matching (e—^ooo, e—^ooe and e—J-oee), 

(3) 

but for all types except e—^ooe values of Xeff 
small and, therefore, these cases are not included into 
Table m 


We paid special attention to calculations of the 
optimal crystal length I and angular detection range 
b^max • 

Note, first, that the total count rate is 
proportional to I in type-I ([^, (151 and independent 
of I in type-II ([^, ( [I^ phase-matching cases. The 
last assertion is correct while in the phase-matching 
inequality (see ([T7| and further) 


|AA:,| = 


aq - 


:79 


27r 


the quadratic term can be omitted. For calcite 
this is true for I ^ Imin = dTry/a^ ^ 0.05 mm . 
But the crystal length can limit the detection angular 
range. For multi-mode detection (see Table 2) we 
assume Omax = 7 r/ 2 ^ [] (we suppose, that all the SPDC 
radiation can be focused on the detector area). But for 
single-mode detection the angular range is defined by 
a gaussian mode divergence dmax = A/ttw, where w is 
the waist of the pump, which should be less than the 
spatial walk-off pl (for the considered crystals p ~ 0 . 1 ). 
So, we obtain Omax o^l/l. It means that in the case of 
type-II phase-matching we need to use as thin crystal 
as possible (we set I = 2lmin = 0.1 mm). In the case 
of type-I phase-matching the crystal length I has to 


be decreased until the integration limits in Eq. (151 


become defined by the angular range. This means that 
the optimal crystal length is 


fa 

/min 

27rc 

UJp 

Up 

27rc 

V/3/ 

. ^min 

3 

3 

^max _ 


Finally, for multi-mode detection and the type-I phase¬ 
matching the crystal has to be taken as long as 
possible (we set 100 mm). 


5. Results and Discussion 

In calculations we took parameters of the most widely 
used single-photon detectors (see Table 2) and the 
most suitable cw lasers: DPSS lasers at 266 and 
532 nm, diode blue-ray laser at 405 nm and HeCd gas 
laser at 325 nm. Typical powers are given in Table B 
It was shown that the use of pulsed pump lasers 
gives no advantages for TOSPDC generation. 

Also we consider a possibility of increasing the 
pump power inside the crystal by means of mirror 
deposition at the rear and front faces of the crystal, 
which turns the latter into a cavity. Accordingly, 
maximal growth of the pump intensity in a crystal 
(cavity) is characterized by the factor e = 1000 .^ 


^ Of course, so big angles are outside of the framework of 
our model but we suppose that, still, it remains reasonable for 
rough estimates . 

® According to |33|, the reflection coefficient of a lossless 
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(3) 

Table 1. Calculated values of effective third-order nonlinear susceptibility Xefp triplets generation rates Rjp, measurement times 
necessary for registration of triple T 3 and double T 2 coincidences for different pump wavelengths Ap and power Wp , type of nonlinear 
media, its length I and type of phase-matching, for different detectors and for the cases with presence (-h) and absence (—) the 
cavity. Easy accessible in an experiment values marked as green, hardly accessible as yellow and unaccessible as red. 


Xp 

(nm) 

Medium 

y(3) 

Xeff 

(lO-i^esu) 

Wp, 

(W) 

Detector 

1 

(mm) 

Cavity 

R'j' 

(Hz) 

Fa 

(days) 

T2 

(days) 

266 

Calcite 
(e—J-ooe) 

0.32 

10 

Si APD 

0.1 

-k 

4.0- 10-3 
4.0- 10-2 

94 

9.4 - 10-2 

15 

1.4- 10-2 

325 

Calcite 
(e—^-ooe) 

0.59 

0.05 

Si APD 
Super Cond. 

0.1 

-k 

1.1 • 10-3 

3.4-10-3 

5 200 

18 000 

750 

1 000 

405 

Calcite 
(e—>ooe) 

0.76 

0.5 

PMD 

Super Cond. 

0.1 

-k 

1 .8- 10--^ 

9.5 • 10-3 

8.1 -10^3 

6 200 

2 .0- 10^1 

370 

532 

Calcite 
(e—J-ooe) 

0.88 

10 

PMD 

Super Cond. 

0.1 

-k 

1 .8- 10--^ 
5.0-10-3 

8.2 -10^3 

1 200 

2 .0- 10^1 
690 

532 

Rutile 
(o—>-eee) 

71.6 

10 

PMD 

Super Cond. 

100 

0.77 

- 

1.5 - 10-2 
8.9-10-^ 

1.1 - 10'^ 
6.6 - 10^ 

2 .8- 10^ 

3.9-103 


Table 2. Characteristics of single-photon detectors to be used for three-photons registration. 


Type 

Xmin XfYiax 
(nm) 

Number 

of spatial modes 

Quantum 
efficiency rj 

Dark count rate Rn^ 
(Hz) 

Jitter St 
(ps) 

Si APD^ 

400-1040 

Multi 

0.1-0.7 

100 

350 

InCaAs APD2 

1000-1650 

Single 

0.1 

3000 

200 

Super Conductive3 

600-1700 

Single 

0.2 

1 

50 

PMlD^^ 

950-1700 

Multi 

0.01 

50000 

70 


1 Excelitas SPCM-AQRH-16 
^ IDQuantique ID210 
3 Scontel SSPD 
Hamamatsu R3809U-69 


Unfortunately, this optimization cannot be used in 
rutile crystals because of their comparatively high 
adsorption. 

In addition to the total triplet generation rates 
Rt, we have calculated and presented in Table 1 the 
time, required for two- and three-photon coincidences, 
( |^ and ([^ at to. 998.00 = 3. . 

One can see, that one of the main problems of 
TOSPDC detection is related to low efficiency and 
high noise of IR detectors. It is really difficult to 
extract a signal from noise even for triple coincidence 
measurements. So we found only one combination 
of parameters when detection of TOSPDC photons 
looks possible. This is the case of a calcite crystal 
with a mirror coverage deposited at rear and front 
faces, the crystal length I = 0.1 mm, and the pump 
parameters Ap=266 nm and IUp=10 W. In this case we 
found T 2 = 20 and T 3 = 135 minutes. 

One more advantage of UV pump is the increase 
of the differential generation rate, because it is 
proportional to u>iOJ 2 U >3 cx Wp (|6 


mirror is Rm Ri 0.999, which gives £ R! (1 — Rm)~^ ~ 1000. 
Losses in mirrors can decrease e making it not higher than 10 
for existing AR coatings. 


Another well-known problem is a low value of 
X^eff For rutile crystal Xgjj is about two orders higher 
if the crystal optical axis is taken parallel to the pump 
propagation direction. But for the phase-matching 
conditions to be satisfied one has to use a crystal with 
periodical poling (quasi-phase-matching). This is an 
evident way for making the type-I and type-II phase¬ 
matching conditions satisfied and TOSPDC photons 

detectable in visible range of wavelengths. Another 

(3) 

possibility of increasing Xeff 1® addition of special 
impurities to crystals which would provide resonance 
enhancement of the third-order susceptibility. 

One more problem is a really multi-mode structure 
of TOPDC generation in bulk crystals, which requires 
multi-mode detection of signals. This problem can be 
solved by producing a wave-guide inside the crystal 
medium, as proposed in |34| . In this way the length of 
interaction can be done arbitrary long. 

6. Conclusion 

Finally, we have performed a detailed analysis of 
three-photon generation in birefringent crystals with 
special attention paid to calcite and rutile crystals. 
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The analysis includes the calculation of differential 
generation rate in unit frequency and transverse wave 
vectors range total count rate ^ for type-I 

and type-II phase-matching and measurement time, 
required for distinguishing signal coincidences from 
noise (( 22 ) and ([2T|). 

The results show that the registration of TOSPDC 
in calcite is possible for the process with the pump at 
266 nm with the presence of a cavity. All the other 
considered cases need too much time for three-photon 
registration. 

This work was supported by the Russian Science 
Foundation (project 14-12-01338). 
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